Scale-dependent rigidity of polymer-ornamented membranes 



(N 
O 

o 

(N 



Laboratoire 



Thomas BickelQ and Carlos M. Marques 
Dynamique des Fluides Complexes, CNRS-ULP, 67084 Strasbourg, France 



We study the fluctuation spectrum of fluid membranes carrying grafted polymers. Contrary 
to usual descriptions, we find that the modifications induced by the polymers cannot be reduced 
to the renormalization of the membrane bending rigidity. Instead we show that the ornamented 
membrane exhibits a scale-dependent elastic modulus that we evaluate. In ornamented lamellar 
stacks, we further show that this leads to a modification of the Caille parameter characterizing the 
power-law singularities of the Bragg peaks. 
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The characterization of polymer-membrane interac- 
tions is a fundamental issue of colloidal science. In cos- 
metics, pharmaceutics or detergency, many formulations 
are suspensions of self-assembled surfactant bilayers with 
polymers added for performance, processing, condition- 
ing or delivery jjj. Likewise, lipid bilayers form the walls 
of living cells and liposomes, and host a great variety of 
macromolecules for coating protection, ion exchange and 
mechanical reinforcement. 0- In many instances, the 
polymers are end-tethered to the soft interfaces. Grafting 
is easily achieved experimentally by using polymer chains 
that carry hydrophobic groups. Typical examples are 
provided by the so-called PEG-lipids, hydrophilic chains 
of polyethylene glycol covalently linked to a double-tail 
phospholipid molecule. Recent studies have shown that 
polymers grafted to bilayers can induce gelation || or 
other phase changes [Q, ^) in liquid lamellar phases. They 
stabilize monodisperse vesicles || and modify the geome- 
try of monolamellar (t], ||] and multilamellar |)| cylindri- 
cal vesicles. They also lead to drastic changes in the 
structure and phase behaviour of ternary amphophilic 
systems [ [Tof . 

As first explained by Canham and Hclfrich, fluid mem- 
branes are fluctuating objects jyj. In thermal equilib- 
rium, the surface assumes all the possible shapes allowed 
by the geometry with the associated Boltzmann proba- 
bility V ~ exp{— Tio/(kBT)}. The Hamiltonian Tio is a 
quadratic function of c\ and C2, the two principal curva- 
tures at any given point of the surface, 



Ho = dS |^(ci + c 2 - 2c ) 2 + kcic 2 | 



(1) 



The spontaneous curvature cq vanishes for symctric bi- 
layers. The Gaussian rigidity R plays an important role 
in the determination of the topology of the system, and 
the integral J dS Rc±C2 does not depend on the particular 
shape of the membrane. The amplitude of the thermal 
fluctuations is controled by the bending rigidity k, that 
ranges from a few k B T to a few tens of k B T . Conversely, 
the analysis of the height correlations in a membrane 
system allows for the determination of the constitutive 
rigidity. Many techniques were specifically developped to 



extract k by comparing experiments and theoretical pre- 
dictions for vesicles, lamellar stacks, bicontinuous phases 
and other geometries. These methods have so far been 
of limited application in polymer-membrane systems be- 
cause of the lack of prediction for the fluctuation spec- 
trum of bilayers in the presence of polymers. In this Let- 
ter, we make a first step to reach this gap by investigating 
the role of grafted chains on the fluctuation spectrum of 
fluid membranes. We first analyze the subtle interplay 
between local monomer concentration and surface curva- 
ture. We show that the chains induce additional height 
correlations that can be interpreted as the result of effec- 
tive, scale-dependent elastic coefficients. For the specific 
case of an ornamented lamellar phase of membranes, we 
also determine the Caille parameter that characterizes 
the power-law singularities around the Bragg peaks. 

In a mean-field description, ornamented membranes 
are assumed to have new parameters K e ff — k + Ak and 
K e ff = R + Ak different from those of the bare mem- 
brane. We consider AT polymers grafted to each side of 
a membrane of area S, the chains being free to diffuse 
along the membrane. The surface coverage o~ = JV/S 
is low enough to avoid polymer/polymer interactions: 
g <C (7* — Hp 2 , with Rp the Flory radius of the chains. 
In this limit of small density, the excess surface energy is 
proportional to a and one anticipates corrections of the 
order of Ak ~ Ak ~ k B T(aRp) in both good and theta 
solvent. The effect of solvent quality and polymer ar- 
chitecture resides in the numerical prefactors of the vari- 
ations. For Gaussian chains tethered by one extremity, 
the calculations can be carried out explicitly and one gets 
Ak = (1 + f ) k B TaR 2 g and Ak = ~2k B TaR 2 g g| ||, 
with R g — (A^a 2 /6) 1 / 2 the radius of gyration of a Gaus- 
sian chain of N monomers of size a. At the overlap con- 
centration a* , the variations are of the order of the ther- 
mal energy k B T. However, for most practical situations, 
the chains can freely move along the membrane surface: 
the degrees of freedom associated with the end positions 
are annealed, as opposed to the quenched results of ref- 
erences |l2|, [l| . A thorough study of the system reveals 
that averaging the free energy over the anchor position is 
not equivalent to averaging the partition function itself. 
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FIG. 1: Diagrammatic representation of the perturbative ex- 
pansion of the Gaussian chain partition function at first order 
(left) and second order (right). 

As a consequence, the response to bending is lowered if 
one correctly accounts for the tendency of the chains to 
be at the top of the membrane ondulation, leading to 

An = k B TaR 2 g , (2) 
Aft = -2k B TaR 2 g . (3) 

The reduction of the effective modulus is due to the linear 
coupling between positions and spontaneous curvature 
and actually follows the general argument first recognized 
by Leibler 

These mean-field results miss however a crucial feature, 
namely that the polymer-membrane interactions strongly 
depend on the length scale. We first consider the parti- 
tion function Zjq of a single Gaussian chain and construct 
a perturbative expansion in the limit of small surface de- 
formations around the flat plane. In absence of overhangs 
and inlets, the surface profile can be described by a single- 
valued function h(r), where r spans the reference plane. 
We expand Z N in a series zffi + + zff + . . ., where 

2$ is of order h l and Z^\z) w z/(^R g ) denotes the 
partition function of a chain tethered at a small distance 
z of a flat, impenetrable wall. It has been shown recently 
that the first-order term of the expansion is related to 
the pressure field applied by the grafted chain to the sur- 
face |l7| . At distance r = \f x 2 + y 2 from the grafting 
point, this entropic pressure decays with the scaling form 
p(r) ~ k B Tr~ 3 (a -C r <C R g ) in both good and theta 
solvents, and then vanishes sharply beyond R g [ fL8| . 

We now evaluate the contribution of the tethered poly- 
mers to the fluctuation spectrum of the membrane. If one 
introduces the Fourier Transform (FT) of the height pro- 
file h(q) — J dre _4qr /i(r), then the energy Ho of the pure 
membrane reads 

H ° = \ J (2^ (q) ^ q) ^ ■ (4) 

Contrary to the previous mean-field results, we assert 
here that the presence of the chains entails non-local 
contributions to the elastic description of the membrane. 
This can be shown by extending the perturbative scheme 
up to second order in h, a calculation that we perform 
analytically for a Gaussian chain. In this case, the single- 



chain partition function obeys the Edwards equation |l9j] 

where (r, z) is the position of the constrained extremity. 
The linearity of eq. (g) implies that each term of the 
development obeys an Edwards equation. The solution of 
the successive orders are coupled through the (Dirichlet) 
boundary condition on the wall Zpj(r, h(r)) = 0. The 
solution of eq. (|^) obeying the impenetrability condition 
is then determined recursively, and is given to each order 
by @ 

2 rN i- /)r>(°) 

Z$(r,z) = YJ Q dn J dr '^p(r,z;^0)z(V,0) , 

(6) 

where the Green function G$ of eq. (pf) satisfies the 
Dirichlet condition on the horizontal plane JTsJ . This so- 
lution corresponds to a "loop expansion" of Z^ '■ the term 
of order i enumerates the conformations that encounter 
the surface i times, as depicted on figure Finally, the 
statistical weight Z^ is averaged for all possible positions 
of the anchoring point along the surface. 

In general, the membrane carries a finite polymer con- 
centration a . We neglect for sake of simplicity the inter- 
actions between the chains, whose origin is twofold: in 
addition to the usual steric repulsions, the mean defor- 
mation resulting from the polymer pressure gives rise to 
attractive, deformation-sharing interactions fl8|| . How- 
ever, the latter do not modify the fluctuation spectrum 
of the membrane. Our results are thus exact in the low- 
density limit crRg <C 1 where direct polymer interactions 
can be neglected. Assuming Af chains on each side of the 
membrane, the total partition function of the system is 
given by the path integral 

2 = m? I m '~ m ° (s / irZ "T 

where the factor l/A/l reflects the fact that particles 
are indistinguishable. The properties of the decorated 
surface are described by the new Hamiltonian 7i e // = 
"Ha + AH, with 

AW = \ J ^U(q)fc(-q)G(q) . (8) 

The algebra will be detailed elsewhere |Q , and we now 
discuss the results. Integrating out the degrees of free- 
dom of the grafted chains generates the additional con- 
tribution 

G(q) = ^{e-^-l + ^} . (9) 
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FIG. 2: Non-local polymer contribution to the fluctuation 
spectrum of a membrane. At the natural length scale R g , the 
function G(q) crosses over from a tension-like behaviour to a 
bending-like behaviour. 



The correlation function G(q) is displayed on figure g. At 
length scales larger than the polymer size qR g <C 1, one 
can extract a correction to the bending rigidity that reads 
An = lim 9 ^ [<7~ 4 G(q)] = k B TaR 2 . We retrieve in this 
limit the mean- field correction (^), as expected. Interest- 
ingly, the polymer-membrane interactions also induce a 
tension-like term at short distances: G(q) ~ 2k B Taq 2 for 
qRg ^> 1. The correlations between chain and membrane 
conformations are stronger inside the polymer mush- 
room: as a consequence, the fluctuations of the surface 
are smoothed out at lenght scales smaller than the poly- 
mer size. The coefficient of this effective surface tension 
corresponds to the two-dimensional perfect gas pressure 
(recall that both sides are ornamented with N chains). 
Note that our results contrast with the usual interfa- 
cial behaviour, where tension is dominant at large scales 
whereas bending effects appear only at shorter scales. We 
should emphasize that the results presented here do not 
depend on the nature of the interface: they can be used 
to study a wide class of surfaces, provided that grafting 
of macromolecules is feasible. 

We now apply our predictions to the lyotropic phase 
L a . This phase consists in a stack of regularly spaced 
membranes, that form a quasi-crystalline structure in one 
direction while retaining their fluid properties in the per- 
pendicular plane. In pure water, the smectic order is 
generally stabilized by long-range electrostatic interac- 
tions. At high salt concentration or for neutral surfac- 
tant molecules, the stability of the phase arises from the 
cntropic repulsions between fluctuating membranes: we 
will refer to this system as to the Helfrich system. In a 
continuum elastic theory, the multilayer membrane sys- 
tem is described by the displacement field u(r) in the z 
direction normal to the layers. The standard theory of a 
pure L a phase is based on the smectic energy density |2l|| 
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(10) 



the subscript j_ referring to the coordinates along the 
layers. K is the smectic curvature modulus and is di- 
rectly related to the membrane bending rigidity through 
the relation K = n/d, with d the average layer spacing. 
The smectic compressibility modulus B reflects the inter- 
lamellar interactions, and is given in a Helfrich system by 
B = 9ir 2 (k B T) 2 /(64Kd 3 ). 

The lyotropic phase is a good candidate to probe 
experimentally the fine modifications to the properties 
of the pure phase induced by the addition of macro- 
molecules [|2j ^3|. X-ray or neutron scattering experi- 
ments allow to measure the exponent r/ that character- 
izes the power-law singularities around the Bragg peak: 
7(g) - \q-q \~ 1+v , with q Q = 2ir/d @. The "strength" 
of the smectic order varies conversely to 77. For pure 
membrane sytems, the FT of the average heigth fluctua- 
tions is (|u(q)| 2 ) = Vk B T/(Bq 2 z + Kq\), which leads to 
the classical prediction first derived by Caille |25| 

, k B T 



(11) 



Likewise, a decorated membrane has according to 
eq. (g) a scale-dependent bending rigidity «(q) = k(1 + 
ag D (q 2 R 2 )) 1 where we define the dimensionless parame- 
ter a = ^°-J%, and g D {x) = £(e~ x -1 + x) is the De- 
bye function p6| . This is certainly the most fundamental 
step for forecasting the properties of membrane-polymer 
complexes, but one should keep in mind that some macro- 
molecular interactions such as repulsion between different 
membranes of the lamellar stack cannot be reduced only 
to the modification of the single membrane parameters. 
Notwithstanding this point, the height fluctuations are 
now given by 

Vk B T 



(Kq)r>, 



Bqt 



Kqi 



(1 + 



ag, 



,(q*R 2 



Kq)l 2 



aKqig D (q 2 ± R 2 ) 



Vk B T 



(Bql 



(12) 



Kql) 

at lowest order in cv. It ensues that the q± corrections 
lead to the same power-law divergence of the scattering 
intensity with new exponent 

k B T 



1 



2 k 



-oR 2 g 



(13) 



The low-g part of G(q) (i.e., the correction to k) actually 
provides the main contribution to the integrals leading 
to (p^). Providing that the interlamellar distance is not 
affected, the tethered chains stiffen the lamellar order, 
that leads to a sharpening of the Bragg peak. A decrease 
of 77 with polymer concentration is indeed observed in 
numerous experiments — see for instance |}| B. How- 
ever, the measured variation of 77 is much sharper than 
the predicted linear behaviour at low coverage. As ex- 
plained above, the chains are also likely to induce non- 
trivial variations of the compressibility B. A simple ex- 
planation of these modifications has been proposed by 
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Castro- Roman et al. H: a reduction in r\ would be due 
to the polymer layer that effectively increases the mem- 
brane thickness. Within this picture, however, it is still 
difficult to explain why the membrane thickness varies 
so sharply with polymer coverage. Our study suggests 
that this could be caused by a non-homogeneous chain 
distribution, arising from the coupling between position 
and curvature. This naturally leads to a higher density 
of polymers in curved regions, which also correspond to 
the contact zones between the fluctuating bilayers. The 
concentration at which rj is significantly modified should 
then be as lower as the coverage needed to roughly have 
one polymer by membrane contact point. 

In conclusion, we have developed a local description of 
the fluctuation spectrum of fluid membranes ornamented 
with grafted polymers. Contrary to current mean-field 
predictions, we find that the effect of the grafted chains 
cannot be reduced to an effective increase of the mem- 
brane bending rigidity. Instead, we show that the elas- 
tic coefficient is a scale-dependent quantity. On scales 
much larger than the polymer size, the bending rigid- 
ity is indeed renormalized. This low-g correction is the 
dominant contribution to the so-called Caille exponent 
measured by X-ray experiments in L a phases. On scales 
smaller than the polymer size, polymer-membrane inter- 
actions smooth the fluctuations, leading to a tension- 
like contribution to the correlation function. Our results 
were obtained for Gaussian chains, a good representation 
for polymers in ^-solvents. It is not clear at this point 
whether or not this description is robust with respect to 
the introduction of monomer-monomer excluded volume 
interactions. In particular, the "tension" regime that de- 
velops at short length scales might exhibit a different, 
non-trivial form. Indeed, the Debye function g D is known 
to scale as g D ~ q~ x ^ , with v = 1/2 in theta solvent 
and v — 3/5 in good solvent conditions. If the relation 
G cx q A g D still holds for chains with excluded volume, 
one should then expect G(q) to behave like G(q) ~ q 7 ^ 3 
for large wavevectors q. Work on this issue is currently 
under progress. 
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